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Abstract In this paper, we discuss self-excited and hidden attractors for systems of 
differential equations. We considered the example of a Lorenz-like system derived from 
the well-known Glukhovsky-Dolghansky and Rabinovich systems, to demonstrate the 
analysis of self-excited and hidden attractors and their characteristics. We applied the 
fishing principle to demonstrate the existence of a homoclinic orbit, proved the dissi- 
pativity and completeness of the system, and found absorbing and positively invariant 
sets. We have shown that this system has a self-excited attractor and a hidden at¬ 
tractor for certain parameters. The upper estimates of the Lyapunov dimension of 
self-excited and hidden attractors were obtained analytically. 


1 Introduction: self-excited and hidden attractors 

When the theories of dynamical systems and oscillations were first developed (see, e.g., the fun¬ 
damental works of Poincare and Lyapunov), researchers mainly focused on analyzing equilibria 
stability and the birth of periodic oscillations . Th e structur es of many applied systems (see, e.g., 
the Rayleigh |148] , Duffing |ffl, van der Pol |144] , Tricomi jl62] , and Beluosov-Zhabotinsky 
systems) are such that it is almost obvious that periodic oscillations exist, because the oscillations 
are excited by an unstable equilibrium. This meant that scientists of that time could compute 
such oscillations (called self-excited oscillations) by constructing a solution using initial data from 
a small neighborhood of the equilibrium, observing how it is attracted, and visualizing the os¬ 
cillation {standard computational procedure). In this procedure, computational methods and the 
engineering notion of a transient process were combined to study oscillations. 

At the end of the 19th century Poincare considered Newtonian dynamics of the three body 
problem, and revealed the possibility of more complicated behaviors of orbits “so tangled that I 
cannot even begin to draw them''. He arrived at the conclusion that ^Ht may happen that small dif¬ 
ferences in the initial positions may lead to enormous differences in the final phenomena". Further 
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analyses and visualizations of such complicated “chaotic” systems became possible in the middle 
of the 20th century after the appearance of powerful computational tools. 

An oscillation can generally be easily numerically localized if the initial data from its open 
neighborhood in the phase space (with the exception of a minor set of points) lead to a long-term 
behavior that approaches the oscillation. From a computational perspective, such an oscillation 
(or set of oscillations) is called an attractor, and its attracting set is called the basin of attraction 
(i.e., a set of initial data for which the trajectories tend to the attractor). 

The first well-known exam ple of a visualization of chaotic behavior in a dynamical system is 
from the work of Lorenz |122 . It corresponds to the excitation of chaotic oscillations from unstable 
equilibria, and could have been found using the standard computational procedure (see Fig. [^. 
Later, various self -excited chaotic attractor s were discovered in many continuous and discrete 
systems (see, e.g.. 


23 26 31 50 123 150 156 





(a) Initial data near the equilib- 
rium So 


(b) Initial data near the equi¬ 
librium Si 


(c) Initial data near the equilib¬ 
rium S2 


Figure 1: Numerical visualization of the classical, self-excited, chaotic attractor in the Lorenz 
system x = 10(1/ — x),y = 28x — y — xz,z = —8/32 -I- xy by the trajectories that start in small 
neighborhoods of the unstable equilibria >S'oq, 2 - Here the separation of the trajectory into transition 
process (green) and approximation of attractor (blue) is rough. 


The study of an autonomous (unperturbed) system typically begins with an analysis of the 
equilibria, which are easily found numerically or analytically. Therefore, from a computationa l 
perspective, it is natural to suggest the following classification of attractors 


80 108 


which is based on the simplicity of finding their basins of attraction in the phase space: 


111 112 


108 111 112^ An attractor is called a self-excited attractor if its basin of at- 


Definition 1 

traction intersects with any open neighborhood of a stationary state (an equilibrium), otherwise it 
is called a hidden attractor. 


The basin of attraction for a hidden attractor is not connected with any equilibrium. For 
example, hidden attractors are attractors in systems with no equilibria or with only one stable 
equilibrium (a special case of the multistability: coexistence of attractors in multistable systems). 
Note that multistability can be inconvenient in various practical applications (see, for example, 
discussions on problems related to the synchronization of coupled multistable systems in |60|[M][W 


142 ). Coexisting self-excited attractors can be found using the standard computational proceduri^ 


^ We have not discussed such possible computational difficulties as caused by Wada and riddled basins 
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whereas there is no standard way of predicting the existence or coexistence of hidden attractors in 
a system. 

Hidden attractors arise in connection with varions fnndamental problems and applied models. 
The problem of analyzing hidden periodic oscillations first arose in the second part of Hilbert’s 
16th problem (1900), which considered the nnmber and mntnal disposition of limit cycles in two- 
dimensional polynomial systems |^. The first nontrivial resnlts were obtained by Bantin (see, 
e.g., [^), which were devoted to the theoretical constrnction of three nested limit cycles aronnd 
one eqnilibrinm in qnadratic systems. Bantin’s method can only be nsed to constrnct nested, small- 
amplitnde limit cycles, which can hardly be visnalized. However, recently an analytical approach 
has been developed, which can be used to effectively visualize nested, normal amplitude limit cycles 


in quadratic systems 75 108 113 




(a) (b) 

Figure 2: Visualization of four limit cycles (green represents stable and red represents unstable) 
in a two-dimensional polynomial quadratic system i = —(aiX^ + bixy + Ciy^ -\- aix -\- /3ij/),y = 
— {a 2 x‘^ + b 2 xy + C 2 y^ +a 2 X + I32y), for the coefficients oi = 6i = /3i = —1, ci = ai = 0, &2 = —2.2, 
and C 2 = —0.7, 02 = 10, a 2 = 72.7778, and /?2 = —0.0015. Localization of three nested limit cycles 
around the stable zero point (green dot) and one limit cycle to the left of the straight line x = —1. 


Later, in the 1950s-1960s, studies of the well-known Markus-Yamabe’s 125 , Aizerman’s |^, 


and Kalman’s conjectures on absolute stability led to the discovery of the possible coexistence 
of a hidden periodic oscillation an d a unique stable stationary point in automatic control systems 
(see [W 15 20 45 79 104 105 143 ; the corresponding discrete examples were considered in |^). 


The Rabinovich system |146J and the Glukhovsky-Dolghansky system |48 

The 


known chaotic systems that have hidden chaotic attractors 172 


interaction of plasma waves and was considered in 1978 by Rabmovich 140 146 


91 


are among the first 
rst one describes the 
Another is a 


model of convective fluid motion and was considered in 1980 by Glukhovsky and Dolghansky 48 
(which we consider in the remainder of this paper). 

Hidden oscillations appear naturally in systems without equilibria, describing various mechan¬ 
ical and electromechanical models with rotation, and elect rical circuits with cylindrical phase 
space. One of the first examples is from a 1902 paper 154 in which Zommerfield analyzed the 





































4 


Will be inserted by the editor 


vibrations caused by a motor driving an unbalanced weight and discovered the so-called Zom- 

42 ). Another well-known chaotic system without equilibria is the 
see also the corresponding Sprott system, which was discovered 


merfield effect (see, e.g 
Nose-Hoover o scillator |53p^ 


independently |156 157 ). In 2001, a hidden chaotic attractor was reported in a power system with 


no equilibria 165 (and references within). 


After the idea of a “ hidden attrac tor” was introduced and the first hidden Chua attractor was 
discovered j^73 76|[M|[T03||111[|112] , hidden attractors have received much attention. Results on 
the study of hidden attractors were presented in a number of invited survey and plenary lectures at 
various international conference^ In 2012, an invited comprehensive surv ey on hidden attractors 


108 


was prepared for the International Journal of Bifurcation and Chaos 

Many researchers are currently studying hidden attractors. Hidden periodic oscillations and 
hidden chaotic attractors have bee n studied in models such as phase-locked loops |68l|7ll, Costas 
loops 16 , drilling systems 110 , DC-DC converters |176 
stabilization systems |^, plasma waves 


175 


aircraft control systems [^, launcher 
interaction , convective fluid motion -EL and many 


others models (see, e.g., (2T|[^[^[^[^[57|[^|6^|6^4| |118}p0|[T^[T34|p^|15^58l|166}]170 


Similar to autonomous systems, when analyzing and visualizing chaotic behaviors of nonau- 
tonomous systems, we can consider the extended phase space and introduce various notions of 
attractors (see, e.g., 25 66 ). Alternatively, we can regard time t as a phase space variable that 


obeys the equation t = 1. For systems that are periodic in time, we can also introduce a cylindrical 
phase space and consider the behavior of trajectories on a Poincare section. 

The consideration of system equilibria and the notions of self-excited and hidden attractors 
are natural for autonomous systems, because their equilibria can be easily found analytically or 
numerically. However, we may use other objects to construct transient processes that lead to 
the discovery of chaotic sets. These objects can be constructed for the considered system or its 
modifications (i.e., instead of analyzing the scenario of the system transitin g int o chaos, we can 
synthesis a new transition sc enari o). For example, we can use perpetual points |l45] or the equilibria 
of the complexified system |l34] . A periodic solution or homoclinic t rajectory can be used in a 
similar way (some examples of theoretical studies can be found in 117|| i27[ 153] ; however the 
presence of chaotic behavior in the considered examples may not imply the existence of a chaotic 
attractor, which can be numerically visualized using the standard computational procedure). 

For nonautonomous systems, depending on the physical problem statement, the notion of self- 
excited and hidden attractors can be introduced with respect to either the stationary states {x{t) = 
xq yt) of the considered nonautonomous system, the stationary points of the system at fixed initial 
time t = to, or the corresponding system without time-varying excitations. If the discrete dynamics 
of the system are considered on a Poincare section, then we can also use stationary or periodic 
points on the section that corresponds to a periodic orbit of the system (the consideration of 
periodic orbits is also natural for discrete systems). 

In the following, we consider an example of a nonautonomous system (a forced Duffing oscilla¬ 
tor), so that we can visualize the chaotic behavior. The classical example of a self-excited chaotic 
attractor (Fig. in a Duffing system x + 0.05i: -I- x^ = 7.5 cos(t) was numerically constructed 
by Ueda in 1961, but it become well-known much later 1163]. To construct a self-excited chaotic 


^ X Int. Workshop on Stability and Oscillations of Nonlinear Control Systems (Russia, 2008), Physics and 
Control |103| (Italy, 2009), 3rd International Conference on Dynamics, Vibration and Control (Hangzhou, 
China, 2010), IFAC 18th World Congress |105] (Italy, 2011), IEEE 5th Int. Workshop on Chaos-Fractals 
Theories and Applications |106] (Dalian, China, 2012), International Conference on Dynamical Systems 
and Applications (Ukraine, 2012), Nostradamus (Czech Republic, 2013) 107 and others. 
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attractor in this system, we use a transient process from the zero equilibrium of the unperturbed 
autonomous system (i.e., without cos(t)) to the attractor (Fig. in the forced system. 



Figure 3: Forced Duffing oscillator: x + O.OSi: + = 7.5cos(<). The {x,x) plane is mapped into 

itself by following the trajectory for time 0 < t < 27r. After a transition process a trajectory from 
the vicinity of the zero stationary point of the unperturbed Duffing oscillator (without 7.5cos(t)) 
visualizes a self-excited chaotic attractor in the forced oscillator. 


Note that if the attracting domain is the whole state space, then the attractor can be visualized 
by any trajectory and the only difference between computations is the timing of the transient 
process. 


2 A Lorenz-like system 

Consider a three-dimensional Lorenz-like system 


X = —a{x — y) — ayz 
y = rx — y — xz 
z = —bz + xy 


( 1 ) 


For a = 0, system Q coincides with the classical Lorenz system 122 . For a > ar and 6=1 after 
a linear change of variables 99 


x —>■ X, y 


C C 

- z, z ^ r - y 


( 2 ) 


system Q takes the following form 


X = —ax + Cz -|- Ayz 
y = Ra - y - xz 
z = —z -I- xy 


( 3 ) 
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with 


C > 0, Ra = 


r{a — ar) 

C 


>0, A=- - ^>0. 

(cr — ar)^ 


( 4 ) 


System (|^ was suggested by Glukhovsky and Dolghansky |^, and describes convective fluid 
motion in an ellipsoidal rotating cavity, which can be interpreted as one of the models of ocean 
flows (see Appendix [A] for a description of this problem). 


In 99 , system (ml was obtained as a linear transformation of the Rabinovich system 146 . It 


describes interactions between waves in plasma |140[ 146]. Additionally, system ([^ describes the 


following physical processes |99]: a rigid body rotation in a resisting medium, the forced motion of 
a gyrostat, a convective motion in harmonically oscillating horizontal fluid layer, and Kolmogorov 
flow. Systems (Q and (|^ are interesting because of the discovery of chaotic attractors in their 
phase spaces. Moreover, system (Q was used to describe the specific mechanism of transition to 
chaos in low-dimensional dynamical systems (gluing bifurcations) |^. 

For system ([^ with cr = ±ar, contains a detailed analysis of equilibria stability and the 
asymptotic behavior of trajectories, and a derivation of the parameter values for which the system 
is integrable. Other researche rs ha ve also considered the analytical and numerical analysis of some 
extensions of system ([l 



_ | |133] . 

Further, following |48J, we consider system ([^ with 

b = 1 , a > 0 , r > 0 , a > ar. 


2.1 Classical scenario of the transition to chaos 


For the Lorenz system 155 , the following classical scenario of transition to chaos is known. Suppose 
that cr and b are fixed (let us consider the classical parameters cr = 10, 6 = 8/3), and that r varies. 
Then, as r increases, the phase space of the Lorenz system is subject to the following sequence of 
bifurcations. For 0 < r < 1, there is globally asymptotically stable zero equilibrium So- For r > 1, 
equilibrium Sq is a saddle, and a pair of symmetric equilibria Si ^2 appears. For 1 < r < r/j « 13.9, 
the separatrices Fi 2 of equilibria Sq are attracted to the equilibria For r = - 13.9, 

the separatrices Ai 2 form two homoclinic trajectories of equilibria S'o (homoclinic butterfly). For 
Th < r < rc ~ 24.06, the separatrices Fi and F 2 tend to S 2 and 5'i, respecti vely. For Tc < r, 
the separatrices F 12 are attracted to a self-excited attractor (see, e.g.. 


155 174 ). For r > 


the 

equilibria 81^2 become unstable. Finally, r = 28 corresponds to the classical self-excited Lorenz 
attractor (see Fig. 0). 

Furthermore, it has been shown that system ([^ follows a similar scenario of transition to chaos. 
However, a substantial distinction of this scenario is the presence of hidden chaotic attractor in 
the phase space of system ([^ for certain parameters values |^ . 

Let us determine the stationary points of system (Q. We can show that for positive parameters, 
if r- < 1, system Q has a unique equilibrium Sq = (0,0,0), which is globally asymptotically 
Lyapunov stable |^. If r > 1, then ([^ possesses three equilibria: a saddle So = (0,0,0) and 
symmetric (with respect to z = 0) equilibria 


Si,2 = i±xi,±yi,zi), 


( 5 ) 


where 


Xi = 


g\/? 

cr -I- a^’ 


2/1 




zi = 


crl -I- ’ 
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and 

a(r — 2) — cr + -\/((T^-arp^P4rnT . 

The characteristic polynomial of the Jacobian matrix of system Q at the point {x, y, z) has the 
form 

xix, y, z) = +pi{x, y, z)X'^ +P 2 {x, y, z)X+P3{x, y, z), 

where 



Pi (a;, y,z)=a + 2, P 2 {x, y, z) = x'^ + ay^ — az^ + (cr + ar)z — rcr + 2cr + 1, 

P 3 {x, y, z) = ax^ + ay^ — az'^ — 2axyz + (ct + ar)xy + (cr + ar)z — ra + cr. 

Following we let ct = 4 and define the stability domain of equilibria S'! 2 - Using the Routh- 
Hurwitz criterion we can obtain the following (see Appendix [b|. 

Proposition 1 The equilibria Si ^2 are stable if 

SaV^ + a(7a - 64:)r^ + (288a + 128)r + 256a - 2048 < 0. (6) 



Figure 4: The stability domain of equilibria 2 for cr = 4. 


The discriminant of the left-hand side of (|6| has only one positive real root, a* « 0.04735. 
So the roots of the polynomial in ^ are as fol ows. For 0 < a < a*, there are three real roots 
ri(a) > r 2 {a) > r 3 (a); for a = a*, there are two real roots: ri(a) and r 2 {a) = r' 3 (a); for a > a*, 
there is one real root ri(a). Thus, for 0 < a < a*, the equilibria Si ^2 are stable for r < r^a) and 
for r 2 (a) < r < ri(a); and for a > a* the equilibria S'i _2 are stable for r < ri(a) (see Fig. ffl. 

Consider the problem of the existence of a homoclinic orbit, which is important in bifurcation 
theory and in scenarios of transition to chaos (see, e.g., [^). For Q and (|^, we can prove the 
existence of ho moclinic trajectories for the zero saddle equilibrium Sq using the fishing principle 
87 96 98 109 . The fishing principle is based on the construction of a special two-dimensional 
manrEoIa such that a separatrix of a saddle point intersects or does not intersect the manifold 
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Figure 5: Separatrices of the saddle S'o = (0,0,0) of system Q for u = 4, a = 0.0052..., and 
r « 7.443. 


for two different system parameter values. Continuity implies the existence of some intermediate 
parameter value for which the separatrix touches the manifold. According to the construction, the 
separatrix must touch a saddle point, so we can numerically localize the birth of a homoclinic orbit. 
A rigorous description is given in Appendix [E| 

For (T = 4, a = 0.0052, and r « 7.443 we numerically obtain a homoclinic trajectory (see Fig.[^. 

We come now to the study of the limit behaviors of trajectories and attractors. We introduce 
some rigorous notions of a dynamical system and attractor and discuss the connection with the 
notions of self-excited and hidden attractors from a computational perspective. 


3 Definitions of attractors 

3.1 Dynamical systems and ordinary differential equations 


Consider an autonomous system of the differential equations 

X = f(x), t S K, x G K", 


(7) 


where f : K" —?> 

The Picard theorem 
that, for any xq G M 


is a continuous vector-function that satisfies a local Lipschitz condition in K". 
(see, e.g., 35 49 ) for a local Lipschitz condition on the function f implies 
, there exists a unique solution x(t, xg) to differential equation (1^ with the 
initial data x(<o,xo) = xq, which is given on a certain finite time interval: t G / C K. Tm theorem 
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35 


49P implies that the solution x(t,xo) 


regarding the continuous dependence on initial data 
continuously depends on xq. 

To study the limit behavior of trajectories and compute the limit values, characterizing trajec¬ 
tories, we consider the solutions of ([^ for t — t -|-cxd or t —> ±00. For arbitrary quadratic systems, 
the existence of solutions for t G [toj +00) does not generally imply the existence of solutions for 
t G (—00, to] (see the classical one-dimensional example i = or multidimensional examples 
from the work on the completeness of quadratic polynomial systems 47 ). It is known that if f is 

(see, e.g., |5^). 


continuously differentiable (f G C^), then f is locally Lipschitz continuous in JK'‘ (see, e.g 
Additionally, if f is locally Lipschitz continuous, then for any xq G K” the solution x(-, xq) : I 
exists on maximal time interval I = (t_, t_|_) G K, where —00 < t_ < u ^ “hOO. If t-j- "hoc, tll6n 
||x(i,xo)|| —)■ 00 for t —)■ and if t- > — oo, then ||x(i,xo)|| —>■ oo for t ^ t- (see, e.g., [TeT] ). This 
implies that a solution of Q is continuous if it remains bounded. For convenience, we introduce a 
set of time values T G The existence and uniqueness of solutions of Q for all t G T can 

be provided, for example, by a global Lipschitz condition. Another effective method for studying 
the boundedness of solutions for alH G T is to construct a Lyapunov function. 

If the existence and uniqueness conditions for all t G T are satisfied, then: 1) the solution of 
([^ satisfies the group property ( |^[4^) 


x(t-I-s,xo) = x(t,x(s,xo)), Vt, sGT, 


( 8 ) 


and 2) x(-, •) : T x K" K" is a continuous mapping according to the theorem of the continuous 
dependence of the solution on the initial data. Thus, if the solutions of ([^ exist and satisfy 
([^ for all t G T, the system generates a dynamical system 17 on the phase space (M”, || • ||). 
Here ||x|| = + ■ ■ ■ + x^ is an Euclidean norm of the vector x = (xi, ..., Xn) G ffi", which 

generates a metric on M". We abbreviate “dynamical system generated by a differential equation” 
to “dynamical system”. Because the initial time is not important for dynamical systems, without 
loss of generality we consider 

x{t,xo) : x(0,a:o) = xq- 

Consider system Q. Its right-hand side is continuously differentiable in K”, which means that 
it is locally Lipschitz continuous i n R" (but not globally Lipschitz continuous). Analogous with the 
results for the Lorenz system 3^ 126], we can prove that the solutions of Q exist for all t G K, 


i.e. system Q is invertible. For this purpose, we can use the Lyapunov function (Appendix]^ 


V{x,y,z) = - 


+ y^ + {a + l)[z- 


(j + r 
fl -f 1, 


> 0 . 


(9) 


Then, system Q generates a dynamical system and we can study its limit behavior and attractors. 


3.2 Classical definitions of attractors 


The notion of an attractor is connected with investigations of the limit behavior of the trajectories 


of dynamical systems. We define attractors as follows 9 19 32 34 82 83 95 160 


^ Similar theorems on the existence, unicity, and continuous dependence on the initial data for solutions 
of system with the discontinuous right-hand side are considered in 64 172 . 
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Definition 2 A set K is said to be positively invariant for a dynamical system if 


and to be invariant if 


x{t, K) CK, Vt > 0, 


x{t, K) =K, Vt > 0, 

where x(t,K) = {x(t,xg) \xq £ K, t> 0}. 

Property 1 Invariant set K is said to be locally attractive for a dynamical system if, for a certain 
e-neighborhood K (e) of set K, 


lim p{K,x{t,xo)) =0, V xq G K{s). 

>- + oo 


Here p{K, x) is a distance from the point x to the set K, defined as 

p{K,x) = inf ||z-x||, 

z^K 


and K{e) is a set of points x for which p{K,x) < e. 

Property 2 Invariant set K is said to be globally attractive for dynamical system if 

lim p{K,x{t,xo)) = 0, V xg € K". 

i—>-+oo 

Property 3 Invariant set K is said to be uniformly locally attractive for a dynamical system if for 
a certain e-neighborhood K{e), any nnmber (5 > 0, and any bonnded set B, there exists a nnmber 
t{6,B) > 0 snch that 

x{t,B n K{e)) C K{S), yt>t{6 ,B). 

Here 

x{t, B n K{e)) = {x(t,xo) I xg G H n K{e)} . 

Property f Invariant set K is said to be uniformly globally attractive for a dynamical system if, for 
any nnmber d > 0 and any bonnded set B C K", there exists a nnmber t{S,B) > 0 such that 

x{t,B)cK{6), yt>t{d,B). 


Definition 3 For a dynamical system, a bounded closed invariant set K is: 

(1) an attractor if it is a locally attractive set (i.e., it satisfies Property^; 

(2) a global attractor if it is a globally attractive set (i.e., it satisfies Property^; 

(3) a B-attractor if it is a uniformly locally attractive set (i.e., it satisfies Property^; or 

(4) o, global B-attractor if it is a uniformly globally attractive set (i.e., it satisfies Property^. 


Remark 1 In the definition of an attractor we assume closeness for the sake of uniqueness. This is 
because the closure of a locally attractive invariant set K is also a locally attractive invariant set (for 
example, consider an attractor with excluded one of the embedded unstable periodic orbits). The 
closeness property is sometimes omitted from the attractor definition (see, e.g., 1^). Additionally, 
the boundedness property is sometimes omitted (see, e.g., j^). For example, a global attractor in 
a system describing a pendulum motion is not bounded in the phase space (but it is bounded 
in the cylindrical phase space). Unbounded attractors are considered for nonautonomous systems 
in the extended phase space. Note that if a dynamical system is defined for t G K, then a locally 
attractive invariant set only contains the whole trajectories, i.e. if xg G K, then x(t,xo) G K for 
Vt G M (see 


32 
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Remark 2 The definition considered here implies that a global B-attractor is also a global attractor. 
Conseqnently, it is rational to introdnce the notion of a minimal global attractor (or minimal 
attractor) |32||34] . This is the smallest bonnded closed invariant set that possesses Property (or 
Property [y. Fnrther, the attractors (global attractors) will be interpreted as minimal attractors 
(minimal global attractors). 


Definition 4 For an attractor K, the basin of attraction is a set B{K) C M" such that 


lim p(itr,x(t,xo)) = 0, 

>-+oo 


V xo G B{K). 


Remark 3 From a compntational perspective, it is not feasible to nnmerically check Property [^for 
all initial states of the phase space of a dynamical system. A natnral generalization of the notion 
of an attractor is consideration of the weaker att ract ion reqnirements: almost everywhere or on a 
set of the positive Lebesgne measnre (see, e.g., |128]). See also trajectory attractors 30 33 151 


To distingnish an artificial compnter generated chaos from a real behavior of the system one can 
consider the shadowing property of the system (see, e.g., the snrvey in |141] ). 

We can typically see an attractor (or global attractor) in nnmerical experiments. The no¬ 
tion of a B-attractor is mostly nsed in the theory of dimensions, where we consider invariant 
sets covered by balls. The nniform attraction reqnirement in Property implies that a global 
B-attractor involves a set of stationary points (5) and the corresponding nnstable manifolds 
kF“(5) = {xo G ffi" I limt_>_oo p(5,x(t, xq)) = 0} (see, e.g., |^|^). The same is trne for B- 
attractor if the considered neighborhood K{e) in Property [^contains some of the stationary points 
from S. From a compntational perspective, nnmerically checking Property|^is also difficnlt. There¬ 
fore if the basin of attraction involves nnstable manifolds of eqnilibria, then compnting the minimal 
attractor and the nnstable manifolds that are attracted to it may be regarded as an approximation 
of B-attractor. For example, consider the visnalization of the classical Lorenz attractor from the 
neighborhood of the zero saddle eqnilibria. Note that a minimal global attractor involves the set S 
and its basin of attraction involves the set kF“(5). Varions analytical-nnmerical met hods for com - 
pnting attractors and their basins of attraction can be fonnd in, for example, ^[39p6|l32[|164Hl77 . 



X 


(a) Two trajectories with symmet¬ 
ric initial data near So 



(b) A trajectory with initial data 
near Si 



(c) A trajectory with initial data 
near S2 


Fignre 6: Self-excited attractor of system ([^ for r = 17, cr = 4, and a = 0.0052, compnted from 
different initial points. 
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4 Self-excited attractor localization 


In 48 system ^ with cr = 4 


was 


studied. Consider the following parameters for system 


0 


cr = 4, a = 0.0052. 


According to Proposition IT if ri « 16.4961242... < r < r 2 ~ 690.6735024, the equilibria S'i _2 
of system Q become (unstable) saddle-focuses. For example, if r = 17, the eigenvalues of the 
equilibria ofsystem Q are the following 

So : 5.8815, -1, -10.8815 

Sia ■■ 0.0084 ± 4.5643i, -6.0168 

and there is a self-excited chaotic attractor in the phase space of system Q. We can easily visualize 
this attractor (Fig.|^ using the standard computational procedure with initial data in the vicinity of 
one of the equilibria 5'o,i,2 on the corresponding unstable manifolds. To improve the approximation 
of the attractor one can consider its neighborhood and compute trajectories from a grid of points 
in this neighborhood. 


5 Hidden attractor localization 

We need a special numerical method to localize the hidden attractor of system ([^, because the 
basin of attraction does not intersect the small neighborhoods of the unstable manifolds of the 
equilibria. One effective method for the numerical localization of hidden attractors is based on a 
homotopy and numerical continuation. We construct a sequence of similar systems such that the 
initial data for numerically computing the oscillating solution (starting oscillation) can be obtained 
analytically for the first (starting) system. For example, it is often possible to consider a starting 
system with a self-excited starting oscillation. Then we numerically track the transformation of 
the starting oscillation while passing between systems. 

In a scenario of transition to chaos in dynamical system there is typically a parameter A G 
[oi, 02 ], the variation of which gives the scenario. We can also artificially introduce the parameter 
A, let it vary in the interval [ 01 , 02 ] (where A = 02 corresponds to the initial system), and choose a 
parameter oi such that we can analytically or computationally find a certain nontrivial attractor 
when A = Oi (often this attractor has a simple form, e.g., periodic). That is, instead of analyzing 
the scenario of a transition into chaos, we can synthesize it. Further, we consider the sequence 
Xj, Ai = oi, Am = 02, Xj G [01,02] such that the distance between Xj and A^+i is sufficiently 
small. Then we numerically investigate changes to the shape of the attractor obtained for Ai = oi. 
If the change in A (from Xj to Aj+i) does not cause a loss of the stability bifurcation of the 
considered attractor, then the attractor for Am = 0,2 (at the end of procedure) is localized. 

Let us construct a line segment on the plane (o, r) that intersects a boundary of the stability 
domain of the equilibria 51^2 (see Fig. 1^. We choose the point Pi(r = 700, a = 0.0052) as the end 
point of the line segment. The eigenvalues for the equilibria of system Q that correspond to the 
parameters Pi are the following: 

So : 50.4741, -1, -55.4741, 

Su2 ■■ -0.1087 ± 10.4543i, -5.7826. 

This means that the equilibria Ai 2 become stable focus-nodes. Now we choose the point Po{r — 
687.5, a = 0.0052) as the initial point of the line segment. This point corresponds to the parameters 
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(a) Po : r = 687.5. Attractor (blue), B-attractor 
(green, blue and red) 


(b) r = 690. Attractor (blue), B-attractor (green, 
blue and red) 



(c) Pi ■. r = 700. B-attractor (green and red) 


Figure 7: B-attractor of system (Q for fixed cr = 4, a = 0.0052, and various r. 


for which system ([^ has a self-excited attractor, which can be computed using the standard 
computational procedure. Then we choose a sufficiently small partition step for the line segment 
and compute a chaotic attractor in the phase space of system (Q at each iteration of the procedure. 
The last computed point at each step is used as the initial point for the computation at the next 
step (the computation time must be sufficiently large). 

In our experiment the length of the path was 2.5 and there were 6 iterations. Here for the 
selected path and partition, we can visualize a hidden attractor of system (Q (see Fig. |^. The 


results of continuation procedure are given in 91 


Note that the choice of path and its partitions in the continuation procedure is not trivial. For 
example, a similar procedure does not lead to a hidden attractor for the following path on the plane 
(a,r). Consider r = 33.51541181, a = 0.04735056... (the rightmost point on the stability domain). 
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Figure 8: Paths [Pq,Pi] and [^ 2 )^ 3 ] in the plane of parameters 
procedure. 


{a, r} used in the continuation 



Figure 9: Hidden attractor (blue) coexist with B-attractor (green outgoing separatrix of the saddle 
S'o attracted to the red equilibria S'i_ 2 ) 


and take a starting point P 2 : r = 33.51541181, a = 0.046 near it (Fig. |^. If we use the partition 
step 0 . 001 , then there are no hidden attractors after crossing the boundary of the stability domain. 
For example, if the end point is P 3 : r = 33.51541181, a = 0.048, there is no chaotic attractor but 
only trivial attractors (the equilibria S*! 2)- 

6 Analytical localization of global attractor via Lyapunov functions 

In the previous sections, we considered the numerical localization of various self-excited and hidden 
attractors of system Q. It is natural to question if these attractors (or the union of attractors) 
are global (in the sense of Definition or if other coexisting attractors can be found. 
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The dissipativity property is important when proving the existence of a bounded global attractor 
for a dynamical system and gives an analytical localization of the global attractor in the phase 
space. The dissipativity of a system, on one hand, proves that there are no trajectories that tend 
to infinity as t —+oo in the phase space and, on the other hand, can be used one to determine 
the boundaries of the domain that all trajectories enter within a finite time. 

Definition 5 A set Bq C M” is said to be absorbing for dynamical system 0 if for any xq S K” 
there exists T = T(xo) such that x(t,xo) S Bq for any t>T. 


Note that the trajectory x(t,xo) with xq G Bq may leave Bq for only a finite time before it returns 
and stays inside for t >T. 


Remark f In 116 the ball Bji = {x S 
if there exists R > 0 such that 


: |x| < i?} was regarded as an absorbing set. In this case. 


limsup |x(t,xo)| < R, 

i—>-oo 


for any xq G 


then it is said that a dynamical system is dissipative in the sense of Levinson . R is called a radius 
of dissipativit'^ 

Definition 6 ( |32 34 ) Dynamical system Q is called (pointwise) dissipativ^ if it possesses a 
bounded absorbing set. 


Theorem 1 ( 

attractor. 


If dynamical system (|^ is dissipative, then it possesses a global B- 


We can effectively prove dissipativity by constructing the Lyapunov function ^|173 . Consider 
a sufficient condition of dissipativity for system ([^. 

Theorem 2 ( [1491171] ) Suppose that there exists continuously differentiable function ^(x) : M" - 
K, possessing the following properties. 

(1) lim|x|_).oo ^(x) = +00, and 

(2) there exist numbers R and k such that for any solution x(t,xo) of system ([^, the condition 
|x(t,xo)| > R implies that y(x(t,xo)) < —>r. 

Then 

(a) any solution x(t,xo) to Q exists at least on [ 0 ,+cx)), so system Q generates a dynamical 
system for any t > 0 and xg € K"; and 

(b) if T] > 0 is such that Bq = {x G K" I C(x) < 77 } D {x G M" | ||x|| < R}, then Bq is a compact 
absorbing set of dynamical system Q. 

More general theorems, connected with the applicat ion of th e Lyapunov functions to the proof 
of dissipativity for dynamical systems can be found in jlOl 


147 


It is known that the Lorenz system is dissipative (it is sufficient to choose the Lyapunov function 
V{x,y,z) = ^{x'^ + y'^ + {z — r — cr) ^)). H owever, for example, one of the Rossler systems is not 
dissipative in the sense of Levinson |115] because the outgoing separatrix is unbounded. In the 
general case, there is an art in the construction of Lyapunov functions which prove dissipativity. 

■* Because any greater radius also satisfies the definition, the minimal R is of interest for the problems 

of attractor localization and definition of ultimate bound. 

5 


Together with the notions of an absorbing set and dissipative system, 


19 


83 


also considered the 


definitions of a B-absorbing set and a B-dissipative system (uniform convergence of trajectories to the 
corresponding B-absorbing set). It is known |19] that if a dynamical system given on (R”, 11 j |) is dissipative, 
then it is also B-dissipative. 
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(a) Absorbing set (grey), self-excited attractor (blue), 
and global B-attractor (blue, green and red) for system 
( |r| with the parameters r — 687.5, cr = 4, and a — 
0.0052. r) ^ 238541.3. 


(b) Absorbing set (grey), hidden attractor (blue), and 
global B-attractor (blue, green and red) for system ^ 
with the parameters r — 700, <7 = 4, and a — 0.0052. 
T) 247230.5. 


Figure 10: Absorbing sets for system Q 


Lemma 1 Dynamical system Q is dissipative. 

The proof is based on Lyapunov function V from ([^ (see Appendix |^ . If i?, ry are chosen as 
in the proof of Theorem Appendix then dynamics system Q has a compact absorbing set 


Bo = 


ix,y,z) : V{x,y,z) 


1 

2 


+ 2/2 + (a + 1) 




For example, for tr > 1, r > 1, a < 1 we can choose R = 


and rj = 2{a + 1)R^ (see Fig. 


10 . 


Note that for system ([^ the ellipsoidal absorbing set Bq can be improved using special additional 
transformations and Yudovich’s theorem (see, e.g., 14 ), similarly to for the Lorenz system. 
There is also a cylindrical positively invariant set for system 0 


C = {\x\<r-\ -r^, y"^ + {z — tY < r"^}, (10) 

a 

because 

( 2/2 + {z — r)2)* < —( 2/2 + {z — r)2) + r2 < 0 Vs, | 2 /| > r or \z\ > 2r 

and 

|s|* < —ctIsI + a\y\ + a| 2 /|| 2 :| <0 |2/| < r, \z\ < 2r, |s| > r + 2 ar^ja. 

Thus, as for the Lorenz system |100], we obtain three different estimates of the attractor: the ball 
Si?, ellipse Bq, and cylinder C. 
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7 Upper estimate of the Lyapunov dimension of attractor 

7.1 Lyapunov exponents and Lyapunov dimension 

Suppose that the right-hand side of system (|^ is sufficiently smooth, and consider a linearized 
system along a solution x(t,xo). We have 


u = J(x(t,xo)) u, u G 


t G 


( 11 ) 


where 


J(x(t,xo)) = 


9/.(x) 

dxj 


lx=x(t,xo) 


is the (n x n) Jacobian matrix evaluated along the trajectory x(t, xg) of system ([^. A fundamental 
matrix X(t,xo) of linearized system 0 is defined by the variational equation 


X{t,xo) = J(x(t,xo))X(t,xo). 


( 12 ) 


We typically set Ar(0,xo) = In, where /„ is the identity matrix. Then u(t,uo) = Ar(t,xo)uo. In 
the general case, u(t, ug) = xo)X“^(0, xgjug. Note that if a solution of nonlinear system (|^ 

is known, then we have 

axg 

Two well-known definitions of Lyapunov exponents are the upper bounds of the exponential 
growth rate of the norms of linearized system solutions (LCEs) |124] and the upper bounds of 
the exponential growth rate of the singular values of fundamental rnatrix of linearized system 

m- 

Let cri(A'(t,Xg)) > • • • > an(X(t,xo)) > 0 denote the singular values of a fundamental matrix 
A'(t,xg) (the square roots of the eigenvalues of the matrix X(t,xo)*X(t,X()) are reordered for each 
t). 

Definition 7 The Lyapunov exponents (LEs) at the point xg are the numbers (or the symbols 
± 00 ) defined by 

LEi(xg) = limsup - ln(Ti(X(t,xg)). (13) 

t—^OO ^ 


(LEs) 


LEs are comm only usecQ in the theories of dynamical systems and attractor dimensions 


19 43 55 85 160 


11 


The LCEs 


124 


and LEs 


>Es [TsT 
|131]. He 


are “often” equal, e.g., for a “typical” system that satisfies the condi- 
However, there are no effective methods for checking Oseledec conditions 


tions of Oseledec theorem 

for a given system: “Oseledec proof is important mathematics, but the method is not helpful in eluci¬ 
dating dynamics” |37| p.118]). For a particular system, LCEs and LEs may be different. For example, 


for the fundamental matrix X{t) 


_ (1 gif) - g f) 


we have the following ordered values: LCEi = 


0 1 

max(limsupA’[y(t)],limsupA’[(ji“^(t)]),LCE 2 = 0; LEi ,2 = max,min(limsupAlgijt)],limsupA’[gr“^(t)]), 

t—>- + oo t—>- + oo t—^ + 00 t —^ + 00 

where X{-) — j log | • |. Thus, in general, the Kaplan-Yorke (Lyapunov) dimensions based on LEs and LCEs 
may be different. Note also that positive largest LCE or LE, computed via the linearization of the system 
along a trajectory, do not necessary imply instability or chaos, because for non-regular linearization there 
are well-known Perron effects of Lyapunov exponent sign reversal |77|[78|[I^. Therefore for the compu¬ 
tation of the Lyapunov dimension of an attractor one has to consider a grid of points on the attractor 
and corresponding local Lyapunov dimensions |81] . More detailed discussion and examples can be found 
in 


74 


102 
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Remark 5 The LEs are independent of the choice o f fnn damental matrix at the point xq unlike 
the Lyapunov characteristic exponents (LCEs, see 124 ). To determine all possible values of LCEs, 
we must consider a normal fundamental matrix. 


We now define a Lyapunov dimension 


62 


Definition 8 A local Lyapunov dimension of a point xq in the phase space of a dynamical system 
is as follows: dim^xg = 0 */LEi(a;o) < 0 and dim^xg = n */X]r=i LE°(a;g) > 0, otherwise 


dim^xg = j(xg) + 


LEi(xg) 


LEj(xg) 


|LEj_|_i(xo)| 


(14) 


where LEi(xo) > ... > LE„(xg) are ordered LEs and j(xg) £ [l,ft] is the smallest natural number 
m such that 


LEi(xg) + ... + LEm(xo) > 0, LEm+i(xg) < 0, 


LEi(xo) + ... + LEm(xg) 


|LE^^i (xg) I 

The Lyapunov dimension of invariant set if of a dynamical system is defined as 

dim^ if = sup dim^xg. 

Xq&K 


< 1 . 


(15) 


Note that, from an applications perspective, an important property of the Lyapunov dimension 


is the chain of inequalities 19 55 56 


dim-r if ^ dim// if ^ dim p K ^ dim p if. 


(16) 


Here dimr if, dim//if, and dim/’if are the topological, Hausdorff, and fractal dimensions of if, 
respectively. 

Along with commonly used numerical methods for estimating and computing the L yapun ov 
dimension, there is an analytical approach that was proposed by Leonov (I^ 94 95 97 99 109 . It 
is based on the direct Lyapunov method and uses Lyapunov-like functions. 

LEs and the Lyapunov dimension are invariant under linear changes of variables (see, e.g., [^). 
Therefore we can apply the linear variable change y = Sx with a nonsingular n x n-matrix S. Then 
system (|^ is transformed into 

Y = Sx = S i{S~^Y) = f(y). 

Consider the linearization along corresponding solution y(t,yg) = Sx{t, S~^xo), that is. 


V = i (y(i, yo)) V, v G K". (17) 

Here the Jacobian matrix is as follows 

J{y{t,yo)) = S J{x{t,xo)) S~^ (18) 

and the corresponding fundamental matrix satisfies E(t,yg) = SX{t,XQ). 

Eor simplicity, let J(x) = Jjxjt, xg) ). Suppose that Ai(x, 5) ^ ^ A„(x, S) are eigenvalues of 

the symmetrized Jacobian matrix (|18|| 


i {SJ{x)S-^ + {SJix)S-^r) . 


(19) 
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Theorem 3 ( | 861|97 ) Given an integer j £ [1,?^] and s € [0,1], suppose that there are a 
continuously differentiable scalar function : K" —> M and a nonsingular matrix S such that 

Ai(x, S') + • • • + Aj(x, S) + sAj_|_i(x, S) + t?(x) <0, \/x £ K. (20) 

Then dim^ K ^ j + s. 

Here d is the derivative of d with respect to the vector field f: 

i?(x) = (grad(t?))*f(x). 

The introduction of the matrix S can be regarded as a change of the space metric. 

Theorem 4 ( [T^[^[^[^ ) A ssume that there are a continuously differentiable scalar function 
d and a nonsingular matrix S such that 

Ai(x,S) + A 2 (x,S) + d(x) < 0, Vxe K". (21) 

Then any solution of system Q bounded on [0, +oo) tends to an equilibrium as t ^ +oo. 

Thus, if ( [2T| holds, then the global attractor of system Q coincides with its stationary set. 
Theorems and give the following results for system Q . 


Theorem 5 Suppose that a > 1. 

If 


(’'+D 


CT\2 ^ 2(cr + 1) 


then any solution of system bounded on [0,+oo) tends to an equilibrium as t ^ + 00 . 

If 






then 


dim^ K ^3 — 


2(a + 2) 


cr + 1 + Y^(cr-l)2+a(f+r)' 


( 22 ) 

(23) 

(24) 


Proof We use the matrix 


S = 


-a-2 0 0 
0 1 0 
0 0 1 


Then the eigenvalues of the corresponding matrix (191 are the following 

A 2 = —1, 


A - ^ ± 1 

“-^ ^ 2 


(cr — 1)^ + a ( 2^ — 


a + ar 


To check property (201 of Theorem and property ( [^ of Theorem we can consider the 
Lyapunov-like function 

^ s)V{x,y,z) 

v(x,y,z) = - -- 

{a-l)‘^ + a{^+ry 
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where 

^ 7 7 9 / o,\ n , , tr + ar 

V(x, y, z) = -x^ + 72 /^ + 7 { 1 + - ) - 2j{r - l)z, 7 = —-— 

(T V cr/ Z(r — Ij 

Finally, for system Q with given S and ■&, if condition ([^ is satisfied and 


s > 


-{a + 3 ) + j(cr - 1)2 + O (f + r)' 


+ l + A/(cr-l)2 + a(f+r)' 


then Theorem 1^ gives (241. If condition ( [^ is valid and s = 0, then the conditions of Theorem]^ 
are satisfied and any solution bounded on[ 0 , + 00 ) tends to an equilibrium as t —> + 00 . □ 

Note that for cr = 4, r = 687.5, and a = 0.0052 the analytical estimate of the Lyapunov 
dimension of the corresponding self-excited attractor is as follows 

diuiLK < 2.890997461... 

and the values of the local Lyapunov dimension at equilibria are 

dimi So = 2.890833450..., dimi, = 2.009763700.... 

Numerically, by an algorithm in Appendix]^ the Lyapunov dimension of the self-excited attractor 
is LD = 2.1405. 

The analytical estimate of the Lyapunov dimension of the hidden attractor for cr = 4, r = 700, 
and a = 0.0052 is as follows 

dimLK < 2.891882349..., 

and the local Lyapunov dimension at the stationary points are the following 

dimLS'o = 2.891767634..., dimi 51,2 = 1.966483617... 

Numerically, the Lyapunov dimension of the hidden attractor is LD = 2.1322. 

Thus, the Lyapunov dimensions of B-attractor (which involve equilibrium Sq) and the global 
attractor are very close to the analytical estimate. 

In the general case the coincidence of the analytical upper estimate with the local Lyapunov 
dimension at a stationary point gives the exact value of the Lyapun ov dimen sion of the global 
attractor (see, e.g., studies of various Lorenz-like systems 90 94 97 99 109 114 ). 
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Appendix 

A Description of the physical problem 

Consider the convection of viscons incompressible flnid motion inside the ellipsoid 


X 3 

— =1, ai > 02 > 03 > 0 

as. 


nnder the condition of stationary inhomogeneons external heating. We assnme that the ellipsoid and 
heat sonrces rotate with constant velocity fJo aronnd the axis. Vector Iq determines the orientation 
of the ellipsoid and has the same direction as the gravity vector g. Vector g is stationary with 
respect to the ellipsoid motion. The valne fJo is assnmed to be snch that the centrifngal forces can 
be neglected when compared with the inflnence of the gravitational field. Consider the case when the 
ellipsoid rotates aronnd the axis 0:3 that has a constant angle a with gravity vector g (|g| = g). The 
vector g is placed in the plane XiXg. Then, SIq = (0, 0, Jig) and Iq = (oi sin a, 0, —03 cos a). Let the 
steady-state temperatnre difference Z\T = (go, 0, 0) be generated along the axis xi (Fig. |ll|. The 
corresponding mathematical model (three-mode model of convection) was obtained by GluEhovsky 
and Dolzhansky 48 in the form (see ([^) 


X = —ax + Cz -f Ayz, 
y = Ra - y - xz, 
z = —z + xy. 


Here 


A 

.4=^1 


Ta =^, Ra = 


A2 


COS a T„ 


C = 


gf^aaQo 

2aia2A/i ’ 

2 afa 2 


—-^5 - T^a sin a, 

03(01 -I- O 2 ) 


, . / /X qPclo cos a , , 

x(!)=^- (“»(*) + 

/,x aPag 


/.X 5/?03 x.x 


and A, /r, /3 are the coefficients of viscosity, heat condnction, and volnme expansion, respectively; 
qi{t), q 2 {t), and 53 (t) ( 93 (t) = 0 ) are temperatnre differences on the principal axes of the ellipsoid; 
a;i(f), uj 2 (t), and 103 ( 1 ) are the projections of the vectors of flnid angnlar velocities on the axes Xi, 
X 2 , and X 3 , respectively. Here 


/,x 5^03 5/303 

The parameters a, Ta, and Ra are the Prandtl, Taylor, and Rayleigh nnmbers, respectively. 


The linear change of variables 48 


X, y^C^y, z^C 
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Figure 11: Illustration of the problem 


transforms system (|^ into the system 


{ i = —ax + z + Acyz, 
y = Rc-y-xz, 
z = -z + xy. 


with 


Rc — RaC, Ac — 


After the linear transformation (see, e.g., 99 ): 


X ^ X, y ^ Rc — 


system (251 takes the form of ([^ with 


RcAc + 1 


z, z ■ 


RcAc +1 


a = 


r, t=^{RcAc + 1). 


{RcAc + iy' a 


(25) 


(26) 


(27) 


B Proof of Proposition 

For system (Q the characteristic polynomial of the Jacobian matrix of the right-hand side at the 
point xg = {xq, 2 / 0 , zq) € has the form 

x(xo) = -bpi(xo)A^ -bp2(xo)A -bp3(xo), 

where 


Pi(xo) =a + 2, 

P 2 (xo) = xl + ayQ — azQ -b (cr + ar)zo — ra + 2a + 1, 

P3(xo) = axl + ayl - azg - 2axoyoZo + {a + ar)xoyo + {a + ar)zo - ra + a. 
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Applying the Hurwitz criterion, the necessary and sufficient stability conditions for stationary 
point xq are the following 


pi(xo) > 0 , 

P 2 (xo) > 0 , 
P 3 (xo) > 0 , and 
Pi(xo)p2(xo) -P3(xo) > 0. 


(28) 

(29) 

(30) 

(31) 


Equilibria S 12 exist for r > 1, and we can check that x(5'i) = x(S' 2 )- For further analysis we 
can introduce 

' (32) 


(33) 


a — arY + dcra — (a — ar)^ > 0. 
ion (291 takes the form 

P 2 ('S'i, 2 ) = {Ci ^/(ct — arY + 4:aa — C'i(cr — ar) — 2a^a{a — ar)'j > 0, 


D = a 

Then for stationary points ([^, condition (291 takes the form 
2 


where 


Cl = a{a — arY + + aa^ > 0 . 

Because a > ar > 0, we have 

Ci\/{a — arY + ‘icra > Ci{a — ar) + 2a^a(a — ar) iff 

2 a^a{a — ar )' ^ 


(cr — ar)'‘ + Aaa > (cr — ar) + 


Cl 


iff 


4:a‘^a{a — ar)"^ Acj‘^a?{a — ar)"^ 

4aa > --C +-- C iff 

Cl > cr(cr — ar)^Ci + a^a{a — ar)^. 

The last inequality is satisfied because 

- (Ci — cr(cr — ar)^Ci — a^a{a — ar)^) = aar{a — ar)^ + a{a^ + ar)^ > 0. 

This implies ( |3^. 

Condition (|30| for 5'i.2 takes the form 

Pz{Si^ 2 ) (a/(ct - ar)'^ + Aaa - (a - ar + 2a)^ ■ 

■ ^(cr — ar)^ + 4(70 — (cr — ar)\J(a — ar)^ + 4cra)^ > 0. 


( 34 ) 


Since 

(cr — ar)^ + 4cra — {a — ar + 2a)^ = 4a^(r — 1) > 0 

and 

(cr — ar)^ + 4cra > (cr — ar)i/Yj^^YrY~-Y4jja iff \/(cr — ar)^ + 4cra > (cr — ar), 

condition ([M) is also satisfied. 
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Condition ( [3T| for stationary points 2 is as follows 

2 


Pi( 5 'i, 2 )p 2 (S'i, 2 ) -P3{Si,2) (c'2\/((t - ar)^ + 4cra- 

—C' 2 (cr — ar) — 2a‘^a{a{a — ar) — 4a)) > 0, 


(35) 


where 


C 2 = ((o'(o’ — ar) — a)^ + a^{ra + 2r — 1) + aa^{a — 2)) . 
If tr > 2, then C 2 > 0 and we can derive a chain of inequalities for ( |3^ : 

C 2 \J[a — ar)^ + 4cra > C' 2 (cr — ar) + 2a^a{a{a — ar) — 4a) iff 


(a — ar)^ + 4cra > I (tr — ar) + 


4cra > 


2a^a{cr{a — ar) — 4a) 
■ ^2 


iff 


4cr^a(cr — ar)(cr((T — ar) — 4a) 4:a‘^a^{a{a — ar) — 4a) 


C2 




iff 


C 2 > cr(cr — ar)(a{<7 — ar) — 4 a)C '2 + cr^a(cr(cr — ar) — 4a)^. 
We can divide the last inequality by (—a^) and rewrite it in the form of polynomial 


x^(T^((T — 2)r^ — a {2a‘^ — Aa^ — 3aa^ + Aaa + 4a) (cr^ + 2(3a — 1)(T^ — 

—8acr + 8a) r — (tr^ + 4a^ — 16a) < 0. 

This inequality corresponds to the stability condition for the equilibria Si 2 - 

□ 


C Proofs of Lemma and the completeness of system (Q 

Suppose that the Lyapunov function has the form 


V{x,y,z) = - 


+y^ + {a + l)[z- 


a -\- r 

d 1. 


(36) 


Here V{x,y,z) —>■ 00 as \{x,y,z)\ —>■ 00. For an arbitrary solution x(t) = (x{t),y(t),z{t)) of 
system Q we have 

V(x, y, z) = x(—ax + cry — ayz) + y{rx — y — xz) + ((a + 1 ) 2 ; — (cr + r)){—z + xy) 

= —ax^ — — (a + l)z^ + (cr + r)z. 

Suppose that e G (0, (a + 1)) and c = min {a, 1, (a + 1) — e} > 0. Then 
V(x, y, z) = —ax^ — j/^ — ((a + 1) — e)z'^ — ez^ + (cr + r)z 


= —ax^ — — {{a + 1 ) — e)z^ — ( Vsz — 


+ , jcr + r) 

2^ J ^ 4e 


„2 , .,2 , 2 n , (cr + r)^ 


< -c(x^ + y + 2: ) + 


46 
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Suppose that > R^. Then a positive >c exists such that 




V(x, y, z) < —cR^ + ' < —>f for > 

4e c 4e 

We choose a number rj > 0 such that 

{{x,y,z) I V{x,y,z) < r?} D {{x,y,z) | < R^} 

i.e., the relation x'^ + y'^ + z"^ < R? implies that 


+ y^ + ia+l){z- 


a + r 
0- “h 1 


= x^ +y^ + z^ + az^ - 2(tT + r)z + 


(a + ry 
n “b 1 


< 277. 


Since 

—2(cr + r)z < 2(cr + r)|z| < 2(cr + r)i?, 
it is sufficient to choose j] > 0 such that 

(a + l)i?^ + 2 (ct + r)i? + < 277 , i.e. 77 >^(a + l)ffl+ ^~*~y ^ . 

a +1 2 \ 0 + 1 / 

Further, we can apply Theorem which implies the Lemma. □ 

Using Lyapunov function p6| , we can prove the boundedness of solutions of system Q for 
t < 0. Note that 


2 z- 


a + r 


+ 


+ lj 2(a + l)2 


(cr + r)^\ / + 3(cr + r)\ 


2(a + 1) / 


2(a + 1) / 


> 0 , 


so the inequality 


2: — 


a + r 


V 

2(a + 1) / 


< 2 2 - 


cr + r \ (o’ + r)^ 

a + 1) 2(a + 1)^ 


is satisfied. This implies that 


U>2a --x2 +2 --7/2 -2(a + l) 


(o’ + ry 


1 


a + r 
a + 1 J ' 4(a + 1)2 
2 


1 (cr + r)2 


= 2a^--x^)+2^--y^)+A^--ia + l)[z- 

1 (cr + 

>2a(-U) + 2(-U)+4(-U)---^ ’ 


a + r 
o + 1 , 


4 a + 1 

1 (cr + ry 
4 a+1 


Suppose that fc = 2cr + 2 + 4, and m = \ 


4 a + 1 

_ 1 (o+rU Tljeji 
4 a+1 ■ 


V + kV > —m. 
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This implies that 

+ ke'^^V > 
at 

Thus for t < 0 we have 

F(0) - e'=V(t) > (me'=* - m)/k 
or 

V{t) < e“^*y(0) + — m)/k. 

This implies that V does not tend to infinity in a finite negative time. Therefore, any solntion 
{x{t),y{t), z{t)) of system Q does not tend to infinity in a finite negative time. Thns, differential 
eqnation Q generates a dynamical system for t S K. 


D Computation of Lyapunov exponents and Lyapunov dimension using 
MATLAB 


The singnlar valne decomposition (SVD) of a fnndamental matrix X(t) has the from 


X(t) = U(t)E(t)V^(t) : = 1 = 


where E(t) = diag{(Ti(t),..., anit)} is a diagonal matrix with positive real diagonal entries known as 
singular values. The singnlar valnes are the sqnare roots of the eigenvalnes of the matrix X(t)*X(t) 
(see [^). Lyapnnov exponents are defined as the npper bonnds of the exponential growth rate of 
the singnlar valnes of the fnndamental matrix of linearized system (see eq. @). 

We now give a MATLAB implementation of the discrete SVD method for compnting Lyapnnov 


exponents based on the prodnct SVD algorithm (see, e.g., 40 159 


Listing 1: productSVD.m - prodnct SVD algorithm 

function [U, R, V] = productSVD(initFactorization, niterations) 

*/, Parameters : 

7« initFactorization - array containing factor matrices of the 
*/, fundamental matrix X, such that: 

7« X = initFactorization(:,:,l) * ... * initFactorization(:,:,end); 

7« niterations - number of iterations in the product SVD algorithm. 

7# dimOde - dimension of the DDEs , nFactors - number of factor matrices 
[", dimOde, nFactors] = size(initFactorization); 

7« A - 2D array of matrices storing the factor matrices at each iteration 
A = zeros(dimOde, dimOde, nFactors, niterations); 

A(:, :, :, 1) = initFactorization; 

7. Q - array of matrices storing orhogonal matrices of the QR decomposition 
Q = zeros(dimOde, dimOde, nFactors+1); 

7# U, V - orthogonal matrices in the SVD decomposition 
U = eye(dimOde); V = eye(dimOde); 

7« R - array of upper triangular factor matrices , such that after 
7. the last iteration \Sigma = R(:,:,l) * ... * R(:,:,end) 

R = zeros(dimOde, dimOde, nFactors); 

7. Main loop 

for ilteration = 1 : niterations 

Q(:, :, nFactors + 1) = eye(dimOde, dimOde); 

for jFactor = nFactors : -1 : 1 
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C = A(:, jFactor, ilteration) * Q(:, jFactor+1); 

[Q(:, jFactor), R(:, jFactor)] = qr(C); 
for kCoord = 1 : dimOde 

if R(kCoord, kCoord, jFactor) < 0 

R(kCoord, jFactor) = -1 + R(kCoord, jFactor); 

Q(:, kCoord, jFactor) = -1 * Q(:, kCoord, jFactor); 

end ; 

end ; 

end ; 

if mod ( i It er at ion , 2) == 1 
U = U + Q(:, :, 1); 

else 

V = V * Q(:, :, 1); 

end 

for jFactor = 1 : nFactors 

A(:, :, jFactor, ilteration + 1) = R(:, :, nFactors-jFactor+1)’; 

end 

end 

end 


Listing 2: computeLEs.m - computation of the Lyapunov exponents 

function LEs = computeLEs(extOde, initPoint, tStep, ... 

nFactors, nSvdIterations, odeSolverOptions) 

7« Parameters : 

*/, extOde - extended ODE system (system of DDEs + var. eq.); 

*/, initPoint - initial point; 

7« tStep - time-step in the factorization procedure; 

7. nFactors - number of factor matrices in the factorization procedure; 

7« nSvdIterations - number of iterations in the product SVD algoritm; 

7. odeSolverOptions - solver options (sover = ode45); 

7# Dimension of the ODE : 
dimOde = length(initPoint); 

7# Dimension of the extended ODE (ODE + Var. Eq.) : 
dimExtOde = dimOde + (dimOde + 1); 

tBegin = 0; tEnd = tStep; 
tSpan = [tBegin, tEnd]; 
initFundMatrix = eye(dimOde); 

initCond = [ initPoint(:); initFundMatrix(:)]; 

X = zeros(dimOde, dimOde, nFactors); 

7« Main loop : factorization of the fundamental matrix 
for iFactor = 1 : nFactors 

[", extOdeSolution] = ode45(extOde , tSpan , initCond, odeSolverOptions); 

X(:, :, iFactor) = reshape(... 

extOdeSolution (end , (dimOde + 1) : dimExtOde), ... 

dimOde, dimOde); 

currInitPoint = extOdeSolution (end , 1 : dimOde); 
currInitFundMatrix = eye(dimOde); 

tBegin = tBegin + tStep; 
tEnd = tEnd + tStep; 
tSpan = [tBegin, tEnd]; 

initCond = [currInitPoint(:); currInitFundMatrix(:)]; 

end 

7. Product SVD of factorization X of the fundamental matrix 
["', R, "1 = productSVD (X, nSvdIterations); 
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*/, Computation of the Lyapunov exponents 
LEs = zerosd, dimOde); 
for jFactor = 1 : nFactors 

LEs = LEs + log (diag(R(:, jFactor))’); 

end ; 

finalTime = tStep * nFactors; 

LEs = LEs / finalTime; 

end 


Listing 3: lyapunovDim.m - computation of the Lyapunov dimension 

function LD = lyapunovDim( LEs ) 

7« For the given array of Lyapunov exponents of a point the function 
7« compute local Lyapunov dimention of this point. 

7« Parameters : 

*/, LEs - array of Lyapunov exponents. 

7« Initialization of the local Lyupunov dimention : 

LD = 0; 

% Number of LCEs : 
nLEs = length(LEs); 

% Sorted LCEs : 

sortedLEs = sort(LEs, ’descend’); 

7. Main loop : 
leSum = sortedLEs (1); 
if ( sortedLEs(l) > 0 ) 
for i = 1 : nLEs-1 

if sortedLEs(i+1) ”= 0 

LD = i + leSum / abs( sortedLEs(i+1) ); 

leSum = leSum + sortedLEs(i+1); 
if leSum < 0 
break ; 

end 

end 

end 

end 

end 


Listing 4: genLorenzSyst.m - generalized Lorenz system Q along with the variational equation 

function OUT = genLorenzSyst(t, x, r, sigma, b, a) 

*/, Generalized Lorenz system with 
"/, parameters: r sigma b a 

OUT(l) = sigma*(x(2) - x(l)) - a*x(2)*x(3): 

0UT(2) = r*x(l) - x(2) - x(l)*x(3); 

0UT(3) = -b*x(3) + x(l)*x(2); 

*/, Jacobian at the point [x(l), x(2), x(3)] 

J = [-sigma, sigma-a*X(3) , -a*x(2); 

r-x(3), -1, -x(l); 

X(2) , X (1) , -b] ; 

X = [x(4), x(7), x(10); 
x(5), x(8), x(ll); 

X(6) , X(9) , X(12)j ; 


"/, Variational equation 
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0UT(4:12) = J*X; 


Listing 5: main.m - computation of the Lyapunov exponents and local Lyapunov dimension for 
the hidden attractor of generalized Lorenz system Q _ 

function main 

*/, Parameters of generalized Lorenz system 
7« that correspond to the hidden attractor 
r = 700; sigma = 4; b = 1; a = 0.0052; 

7« Initial point for trajectory that visualizes the hidden attractor 
xO = [-14.551336132013954 -173.86811769236883 718.92035664071227]; 

tStep = 0.1; 
nFactors = 10000; 
nSvdIterations = 3; 

7« ODE solver parameters 

acc = le-8; RelTol = acc; AbsTol = acc; InitialStep = acc/10; 
odeSolverOptions = odeset(’RelTol’, RelTol, ’AbsTol’, AbsTol, ... 

’InitialStep’, InitialStep, ’NormControl’, ’on’); 

LEs = computeLEs(0(t, x) genLorenzSyst(t, x, r, sigma, b, a), ... 

xO, tStep , nFactors, nSvdIterations, odeSolverOptions); 

f printf (’Lyapunov exponents: 7,6.4f, 7#6.4f, 7#6.4f\n’, LEs); 

LD = lyapunovDim(LEs); 

f printf (’Lyapunov dimension: 7#6.4f\n’, LD); 
end 


E Fishing principle and the existence of a homoclinic orbit in the 
Glukhovsky-Dolghansky system 

E.l Fishing Principle 

Consider autonomous system of differential equations Q with the parameter 

X = f(x, q), t G R, X S R", q G R"*. 


(37) 


Let 7 (s), s G [0,1] be a smooth path in the space of the parameter {q} = R"*. Consider the following 
Tricomi problem |162 : Is there a point q^ G j(s) for which system ( |37[ ) with go has a homoclinic 
trajectory? _ 

Consider system (37l with g = 7 ( 5 ), and introduce the following notions. Let x(t, s)+ be an 
outgoing separatrix of the saddle point xq (i.e. lim x(t, s)+ = xq) with a one-dimensional unstable 

i—>• —00 

manifold. Define by xr 2 (s)+ the point of the first crossing of separatrix x(t, s)"*" with the closed set 
12 : 

x(t, s)"*" Gi2, t G (— 00 , T), 
x(r, s)+ = xr 2 (s)+ G 12. 

If there is no such crossing, we assume that xr 3 (s)+ = 0 (the empty set). 
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Theorem 6 (Fishing Principle |87 96 98 ) Suppose that for the path 7 ( 5 ) there is an (n — 1)- 
dimensional bounded manifold f2 with a piecewise-smooth edge dfl that possesses the following 
properties. 

1. For any x G f2\ dfl and s G [0,1], the vector /(x, 7 ( 5 )) is transversal to the manifold 17 \ 917, 

2. for any s G [0,1], /(xo, 7 (s)) = 0, the point Xq G 917 is a sa ddle; 

3. for s = 0 the inclusion Xf 2 ( 0 )^ G 17 \ 917 is valid (Fig. 12a), 

4 . for s = 1 the relation Xi 7 (l)^ = is valid (i.e. Xf 2 (l)+ is an empty set), 

5. for any s G [0,1] and y G dfl \ Xq there exists a neighborhood U{y,d) = {x| |x — y| < (5} such 

that xa{s)~^GU{y,6). 

// cond*9ons[^-l^ are satisfied, t hen t here exists sq G [0,1] such that x(t, Sq)'^ is a homoclinic 
trajectory of the saddle point Xq (Fig. 12b). 


x(t,sf 




[0,so). 


Figure 12: Separatrix x(l, s) + , where s G [0,so]. 


The fishing principle can be interpreted as follows. Fig. |12a| shows a fisherman at the point xq 
with a fishing rod x{t, s)'*'. The manifold 17 is a lake surface and 917 is a shore line. When s = 0, a 
fish has been caught by the fishing rod. Then, x(t, s)+, s G [0, sq) is a path taken by the fishing rod 
when it brings the fish to the shore. Assumption 5) implies that the fish cannot be taken to the 
shore 917\xo, because 9l7\xo is a forbidden zone. Therefore, only the situation shown in Fig. 12b 
is possible (i.e., at s = sq the fisherman has caught a fish). This corresponds to a homoclinic orbit. 

Now let us describe the numerical procedure for defining the point F on the path 7 ( 5 ), which 
corresponds to a homoclinic trajectory. Here we assume that conditionsQ,]^, andl^ of the fishing 
Principle are satisfied. Consider a sequence of paths 7 j(s) C {7j_i(s)7 s G [0,1]} C { 7 ( 3 ), s G 
[0,1]}, Vs G [0,1] such that the length { 7 ^( 3 )} tends to zero as j -G +(X). Condition 3) is satisfied 
for 7 j( 0 ) and condition 4) is satisfied for 7 ^ (1). This sequence can be obtained if the paths 7 and 
7 j are sequentially divided into two paths of the same length and we choose the path such that for 
its end points condition 3) is satisfied and condition 4) is not satisfied (or vice versa). Obviously, 
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the sequence 7 j(s), s G [0,1] is contracted to the point r G { 7 j(s), s G [0,1]}, Vj. This point 
corresponds to a homoclinic trajectory of system (371. 

Now, consider the conditions of the non-existence of a homoclinic orbit. Consider the Jacobian 
matrix of system ( [3^ 

= ^(x,7(s)). 

Let Ai(x, s. S') ^ ^ A„(x, s, S) denote the eigenvalues of the symmetrized matrix 


-(SJ(x,s)S-i + (SJ(x,s)S-i)*), 
where S is a nonsing ular matrix. 

Suppose system (371 has a saddle point xq = xo(s). Vs G [0,1], the point xq belongs to a 
positively invariant bounded set K, and J{xq,s) has only real eigenvalues. 

Theorem 7 ( dzi) Assume that there are a continuously dijferentiable scalar function V(x, s) 
and a nonsingular matrix S such that for system (37l with q = 7 ( 3 ), the inequality 

Ai(x, S)-I-A 2 (x, S) -I-ii(x) < 0, Vx G iV, Vs G [0,1] (38) 

is satisfied. Then system ( |37| has no homoclinic trajectories for all s G [0,1] such that 

lim x(t) = lim x(t) = xq. (39) 

>• —00 t—>-+oo 


E.2 Existence of a homoclinic trajectories in the Glukhovsky-Dolzhansky system 

Consider the separatrix x'^(t), y'^(t), z~^{t) of the zero saddle point of system (Q, where x{f)~^ > 0, 
yt G (— oo,t), r is a number, and lim x(t)~'~ = 0 (i.e. positive outgoing separatrix is considered). 

t—>• —00 

Define the manifold f2 as 

f2 = {x = 0, y < 0, y{a — az) < 0, -I- < 2r^}. 

Check condition 1). 

Inside the set Q\dfl we have 

X = y{a — az) < 0. 

Check condition 5). 

a) On Bi = {x = 0, y = 0, —y/^r < z < cr/a} system (Q has the solution 

x{t) = y{t) = 0, z{t) = z(0) exp(-t). 

b) On B 2 = {x = 0, y < 0, z = a/a, y^ + z^ < 2r^} we have 


X = ay. 

Therefore the local behavior of trajectories in the neighborhood of B 2 is shown in Fig. 

c) The set B 3 = {x = 0, y < 0, —\/2r < z < a/a, y"^ + z^ = 2r^} is located outside of the 
positively invariant cylinder C (see eq. (H)) . Thus, the separatrices of the zero saddle point (which 
belongs to C) can not reach the set B^. 
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Figure 13: Manifold 1?. 



Figure 14: Local behavior of the trajectories of system Q in the neighborhood of set B 2 (cr = 4, 
a = 0.52, r = 10'^) 


Check condition 3). 

Consider the development of the asymptotic integration of system Q |^. Assume that 

ar = c — Xs + O(s^), (40) 

where c and A are some numbers and £ = l/y^ is a small parameter. 

Lemma 2 For any a > c, a > 1 there exists a time T > 0 such that for sufficiently large r, 
{T)^y^ {T)^ z'^ {T)) € fl\df2 (i.e. condition^^ of the fishing principle is valid). 

Proof (sketch). Using the transformation 

t^—j=, X ^ \/rx, y ^ ry, z ^ rz, (41) 

y/r 

we can obtain 

X = ay — eax — (c — Ae + 0{e'^))yz 

y = X — ey — xz (42) 

z = —ez + xy. 
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a) Consider the zero approximation of system ( [4^ (system ( [4^ withont Xe and e = 0) and its 
solntion {xo{t),yQ{t), zo{t)). There are two independent integrals 

V{xo{t),yo{t),zo{t)) = {a - c)zo{tf + ayoit)'^ - Xo{tf = Ci, 

W = yo{tf + zo(i)^ - 2zo(t) = C' 2 . 


Thns, the positive ontgoing separatrix (t), j/,| (t), (t) of the saddle point {x = y = z = Q) oi 

zero approximation of system (421 belongs to the intersection of snrfaces V = 0 and VC = 0, i.e. 

V{x^it)>yoit),Zo{t)) = 0 = W{x^{t),y+{t),z^{t)),yt G {- 00 ,+ 00 ). (43) 


From (43 1 it follows that Xq (t) 


^ 0,yt G K" and 


lim 

>- + 00 


x^{t) 


lim 

t —f-|-oo 


vUt) = 


lim z^ (t) = 0. 

(->■+00 ^ 




(b) First approximation of system (42 1 : for 
xt (^) ~ 0 the positive outgoing separatrix is 
located inside magenta domain i?, and near red 
circuit and the border of blue ellipsoid (48 1 


Fignre 15: cr = 2.3445, a = 0.0065, r = 300, c = 2, A = 1. 


b) Consider the first approximation of system ( [3^ (system ( [42| withont 0(e^)). For the small 
valnes of e the ontgoing separatrix (x^ {t),y^ (t), of the zero saddle point of the first approx¬ 

imation of system ( [42| is close to {xq (t), yQ (t), Zq (t)) on (— oo,t). Therefore for snfficiently small 
e and some r the separatrix {xf (t), y^ (t), z^ (t)) reaches 5-vicinity of the zero saddle. Then there 
exists finite r = r(£, 6 ) snch that 

|xf (r(e,(5))| < <5, \yt{T{e,S))\ < 6, |z)^(r(£,(5))| < <5. 

Consider two fnnctions 


14 (x, y,z) = {a - c + Xe)z‘^ + ay"^ - x^ 


W = y^ + z^ - 2z. 


(44) 
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For the derivatives of ([44]| along the positive outgoing separatrix we have 


W. Vi W) = -‘^sVeixf {t),yt{t),zt (t)) + 2£(cr - l)x"[{tf, 
^W{t) = ^W{xt{t),ytit),zt{t)) = -2eW{xt{t),yt{t),z+{t))-2ezt{t). 


( 45 ) 


dt 

Integrating 

we obtain 


dt 

from —oo to t and taking into account 

lim I4(ro) = lim VF(ro)=0, 

To —>• — OO To —>• — OO 


r 

Vs{t) = 2 e{a — 1) j dt = 2e{a — 1 )Mq +o{e), 

— OO 

r 

kF(T) = — 2£ J (s))'^ ds =—2eNo + o{e), 

— OO 

where Mq and Nq are some positive numbers. If zi and Z 2 satisfy 

y,(O,O,0i) = 2£(a- l)Mo, IF(0,0 ,z2) = -2eNo, 

then 

l2Mo{a - 1) _ _ 2No 


Zi = 


a — c + Xe 


Z2=e 


l + y'l-2Noe' 


(46) 

(47) 

(48) 

(49) 


Hence the situation shown in Fig. 15b occurs in the neighborhood of the saddle point {x = y = 
z = 0 ) for the surfaces 


Ve{x,y,z) = 2s(a- l)Mo, W(x,y,z) = 2sNo. 


(50) 


The separatrix (x^(T),y^(T), z^(t) has to be near the surfaces Ve{x,y,z) and W{x,y,z). From 
the mutual disposition of surfaces (501 and different order of smallness in ( [4^ it follows that if 
xfir) > 0, then yf (t) < 0 and xf (r) < 0 for sufficiently small £. Moreover, if (t) < 0 for t > t 
and xf (t) < 0. This implies the existence of T > r such that xf (T) = 0. We can obtain similar 
results for the case a;^(r) < 0 (then j/)~(t) > 0). The behavior of separatrix {x~^ (t) , y~^ (t) , z^ (t) 
is consistent with the behavior of (a:^(r), yf (t), z f ( r). 

Check condition 4). We now check condition^ for system (Q with parameters ^7\ . For this 
system it was proved (T^ pp. 276-277, 269-272] that if 


Rc < 


4(cr + 1) 


(51) 


1 + y/l + SAci^cr + 1) 

then condition (381 of Theorem ^is satisfied for S = I, V{x, s) = 0. 

Now we can show that if (|M| holds, then condition for system ([^ is also satisfied. Consider 
the path 

Reis), Ac(s), (t(s) = a, 

A,(0) = 0, 71,(1) = 71„ 

r, / \ ( 2ct 4((T + 1) \ 

y 1 + y/l + 4Acis)a 1 + -y/l + 871 c(s)(ct + 1) J 
Rei0) = a{l + S), 


( 52 ) 
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89 96 ). If for some si S [0,1] condition 


where <5 is a small positive nnmber. 

For s = 0 condition |4]| is satisfied (see, e.g. 
is not satisfied, then condition ^ is satisfied for Si- In this case Theorem implies that there 
exists S 2 G [0, si] for which a homoclinic trajectory exists. Bnt i?c(s) is chosen in snch a way 
that conditions of Theorem]^ are valid and hence the homoclinic trajectories do not exist. This 
contradiction proves the fnlfillment of condition]^ of Theorem]^ for all s € [0,1]. Conditionis 
checked. 

Check condition 2). From (|5^ it is obvions that condition 2 is satisfied. 


Remark 6 For c = a Lemma is not valid since a positive ontgoing separatrix of the zero approx¬ 
imation of system (42 I follows a heteroclitic orbit 


lim j;[](t) = lim j/°(t) = lim z°(t) = 2. 

t —^~|”00 t —^~1“00 t —^~1“C50 


In this case we may consider a seqnence of systems close to ([^. For example, instead of 


consider a = a(/3fc) = Jly , where /?*, are a small positive nnmbers and lim /3fc = 0, snch 
that path ([5^ satisfies condition 4) of the fishing principle. Then, nsing Lemma and the fishing 


principle, we get the seqnences of and corresponding homoclinic orbits. Choosing a convergent 


snbseqnence from and nsing Arzela-Ascoli theorem, we can jnstify the existence of a homoclinic 
orbit in the initial system with a(0). 

Note also that since a and r are varying in the asymptotic integration, dH is also varying. 


Finally we get the following 

Theorem 8 For any fixed > 0, cr > 1 there exists a number 

^ f 2a \ 

Rc G [ : ~ ’ , “boo I 

V1 T Vl -b 4 Ac(t J 

such that system Q with parameters (|27|l has a homoclinic trajectory of the zero saddle point. 
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